Integral representations for nonperturbative generalized parton distributions in terms 

of perturbative diagrams 
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An integral representation is suggested for generalized parton distributions which automatically 
satisfies the polynomiality and positivity constraints. This representation has the form of an integral 
of perturbative triangle diagrams over the masses of three propagators with an appropriate weight 
depending on these masses. An arbitrary D term can be added. 
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I. INTRODUCTION 

Generalized parton distributions (GPDs) [1|, 
IE S 0- EJ play an important role in the QCD anal- 
ysis of various hard phenomena such as deeply virtual 
Compton scattering and hard exclusive meson produc- 
tion. GPDs are defined in terms of nondiagonal hadron 
matrix elements of two quark (gluon) fields separated by 
a light-like interval. GPDs contain a vast amount of non- 
perturbative information about the quark-gluon struc- 
ture of hadrons. In particular, the usual forward parton 
distributions (FPDs) and the form factors of hadrons can 
be expressed via GPDs. 

In contrast with form factors and FPDs which can 
be directly accessed experimentally, the case of GPDs 
is much more involved: the experimental data can pro- 
vide information only about some integrals containing 
GPDs. On the theoretical side, GPDs are typical non- 
perturbative quantities. Although there are no reliable 
methods for the calculation of GPDs from the first princi- 
ples of QCD, still the theory imposes certain constraints 
on GPDs which should be taken into account in the anal- 
ysis of the experimental data. Among the general theo- 
retical constraints on GPDs an important role is played 
by the polynomiality of the Mellin moments 01 and by 

thepositivity bounds E3, E3, E3, EM EH E3, EM Ea 
mEllilEI. 

Since the GPDs cannot be measured directly, in the 
analysis of experimental data one has to deal with mod- 
els of GPDs. It would be preferable to use those models 
which are compatible with the polynomiality and posi- 
tivity constraints. However, the realization of this idea 
meets with problems. The polynomiality property holds 
automatically if one uses the so-called double distribution 
representation for the GPDs HUH- 

However, the dou- 
ble distribution representation does not guarantee posi- 
tivity. The general solution of the positivity bounds is 
also known |24j but the problem remains how to describe 
the class of models of GPDs which satisfy both positivity 
and polynomiality conditions. 

In Ref. [24[ an ansatz for GPDs is suggested which 



automatically obeys both positivity and polynomiality 
constraints. The method of Ref. |2J] is based on a for- 
mal mathematical construction rather than on physical 
arguments. In this paper another approach is taken. We 
start from an analysis of simple perturbative graphs for 
GPDs. On general grounds these graphs must obey both 
positivity and polynomiality constraints. We check the 
positivity explicitly. Next we notice that the set of func- 
tions obeying both polynomiality and positivity condi- 
tions is convex. Therefore taking linear combinations of 
perturbative graphs for different theories weighted with 
positive coefficients we obtain new solutions of the pos- 
itivity and polynomiality constraints. The words "dif- 
ferent theories" mean that we can average over various 
parameters: masses, vertices, couplings, sets of fields, etc. 
At first sight this approach looks like an artificial trick 
rather than physics. However, in this paper we reveal 
certain structures standing behind the leading-order per- 
turbative graphs for GPDs in various theories and show 
that these structures can be used as a sort of elemen- 
tary blocks for the construction of a rather wide class of 
models of GPDs obeying both polynomiality and positiv- 
ity constraints. The analysis is restricted to the case of 
spinless hadrons (e.g. pions) but the methods suggested 
here allow a straightforward generalization for the more 
interesting case of nucleon. 



The structure of the paper is as follows. Sec. [n] de- 
scribes notations used for GPDs in the usual and im- 
pact parameter representations. Sec. IIIII contains a brief 
review of the positivity and polynomiality properties of 
GPDs. Sec. IIVI explains how perturbative GPDs can be 
used in order to construct solutions of the polynomiality 
and positivity constraints on GPDs. In Sectionslyland lvTl 
the leading order perturbative GPDs are analyzed in the 
cj) 3 and Yukawa models respectively. In Sec. IVIl| integral 
representations are suggested for GPDs which automat- 
ically obey the polynomiality and positivity constraints. 
In Sec. IVIIII the consistency of the approach is tested 
by checking the positivity of the corresponding forward 
parton distributions. 
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TABLE I: Light-ray operators 0^ N ' (X,n) for various types of 
partons and the corresponding parameter N. 



II. GENERALIZED PARTON DISTRIBUTIONS 

The GPDs#W(x,f,t) 

can be denned in terms of ma- 
trix elements of parton light-ray operators fjW over the 
hadron states P&) with momenta P&: 



dX 
2tt 



exp 



(iAx)(P 2 |C»W(A,n)|P 1 ). (1) 



The light-like vector n, 



= 0. 



(2) 



is normalized by the condition 



n(Pi+P 2 ) = 2. (3) 

We use the standard notation of Ji |l4j for parameters 
A, t and £: 

1 

2 



A = P 2 -P 1 , £ = -i(nA), i = A 2 



(4) 



The definitions of the light-ray operators (A, n) for 
various types of partons are listed in Table [I] We have 
included the scalar field 4> in this table since the GPD 
of the <fi 3 model will be an essential ingredient of our 
construction. The last column of this table contains the 
number N of factors n M appearing in the light-ray op- 
erator 0^ N ' (A,n). This number N plays an important 
role in the formulation of the positivity bounds and of 
the polynomiality conditions and we include N in the 
notation (JTJ of the GPD H^(x, f, t). 

In the frame where (Pi + Pa) = and n =0, 
the transverse component A^ of the hadron momentum 
transfer A Q is connected with the parameter t = A 2 
by the following relation 



t = -- 



|A J 



4£ 2 M 2 



1 



(5) 



Below for the analysis of the positivity bounds we shall 
need the impact parameter representation for GPDs 0, 
13 El 111 HH|. We define the GPD in the impact 
parameter representation as follows: 



p (JV) M,& J 



ri 2 A J 



IA 



±12 



exp ^(A^)] 



4£ 2 Af 5 



(6) 



III. POLYNOMIALITY AND POSITIVITY 

In this section we briefly describe the polynomiality 
and positivity properties of GPDs. The polynomiality 
means that Mellin moments in x of GPD HW(x,£,t), 



dxx m H<- N \x^,t) =P m+N (£,t), 



(7) 



must be polynomials in £ of degree m + N . 

The positivity bounds on GPDs have a simple form in 
the impact parameter representation @. In Refs. [221I23T ] 
the following inequality was derived: 



d£ / dx(l - x) 



1-x 
~ 1 



1-x 



x F 



(AO 



i-e 



> 0. 



(8) 



This inequality was derived in Refs. |22J, |23j for the case 
N — 1 and the generalization for arbitrary N is straight- 
forward. 

The inequality JSJl should hold for any function p(z). 
Therefore we deal with an infinite set of positivity bounds 
on GPDs. The inequality (JHJ) (with its generalizations 
for the nonzero-spin hadrons and for the full set of the 
twist-two light-ray operators) covers various ine qual ities 
suggested for GPDs [H [H Q G3 El [H M, M, 
I20I |21| as particular cases corresponding to some special 
choice of the function p{z). 

It is well known that the double distribution represen- 
tation HHI3 with the D term H3 



H(x,i,t)= J dad(35{x 

M + I/3|<1 



■Za-0)F D (a,0,t) 



+ 0m-\x\)D(-,t)sign(O 



(9) 



guarantees the polynomiality property Q. Another in- 
teresting parametrization for GPDs supporting the poly- 
nomiality was suggested in Ref. [28j . 

On the other hand, as shown in Ref. [24| (see also Ap- 
pendix [SJ, the positivity bound on GPDs © is equiva- 
lent to the following representation for GPDs in the im- 
pact parameter representation in the region x > |£|: 



N+l 



Eo»(i+f' (1 - 0frJ 



1-x 

~ 1 



( 10 ) 

with arbitrary functions Q n . Instead of the discrete sum- 
mation over n one can use the integration over continuous 
parameters. 
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Although both polynomiality and positivity are basic 
properties that must hold in any reasonable model of 
GPDs, usually the model building community meets a 
dilemma: one can use the double distribution represen- 
tation @ but it does not guarantee that the infinite set of 
inequalities JBJ will be satisfied [2fj liffll l3ll | . Alternatively 
one can build the models based on the representation l|l(J|) 
or on the so called overlap representation |15| , which also 
automatically obeys the positivity bounds, but then one 
meets problems with the polynomiality. In this paper 
a rather general representation for GPDs is suggested 
which guarantees both positivity and polynomiality. 



IV. GENERAL METHOD 

One could consider the construction of a representa- 
tion for GPDs which solves simultaneously positivity and 
polynomiality constraints as a pure mathematical prob- 
lem, looking for functions Q n <|K>P which allow the double 
distribution representation ©: 



(i + e)fc J 



= /|^ r exp[i(A^)] J dadp 

\a\+\/3\<l 



x S (x - - £a) F% (a,P 



(ii) 



If we manage to find a large set of such functions Q„, 
then taking linear combinations with positive coefficients 
we can construct many solutions of the positivity and 
polynomiality constraints. This strategy was used in Ref. 

On the other hand, the solution of the positivity and 
polynomiality constraints is a physical problem and in- 
stead of using formal mathematical methods one can try 
to solve this problem relying on physical arguments. The 
polynomiality and positivity constraints hold in any rea- 
sonable quantum field theory. In particular, we expect 
these properties in the leading order perturbative dia- 
grams for GPDs in various field theories. Now it makes 
sense to notice that the form of the polynomiality and 
positivity constraints is sensitive to the spins of par- 
tons and hadrons but not to the dynamics of the theory. 
Therefore taking a formal "superposition" of the leading- 
order perturbative GPDs Hm (x, £, t) over various models 
Ai (and over various values for the parameters of these 
models) with arbitrary positive coefficients cvi, 

H{x,£,t) =J2 c MH M {x,Z,t), cm>0, (12) 

M 

we also obtain a representation for GPDs which auto- 
matically obeys both polynomiality and positivity con- 
straints. 



At first sight the mathematical approach based on re- 
lations 110|l and the diagrammatic method (|12p are abso- 
lutely different ways to solve the positivity and polynomi- 
ality constraints. But there is a deep relation between the 
two approaches. The leading order perturbative GPDs 
HM{x,£,t) obey the positivity condition JSJ. Therefore 
these perturbative GPDs Hm can be represented in the 
form HI U| l in the impact parameter representation JfJJl. 
Actually the decomposition IjlOp arises automatically if 
one computes the leading order triangle Feynman dia- 
grams for Fm ( x, 6 , b- 1 ) directly in the impact parameter 
representation 32] . The sum over n on the right-hand 
side (RHS) of Eq. ltTU|) is nothing else but the sum over 
the types and polarizations of the intermediate particles 
of the triangle graphs for GPDs. Therefore this sum is 
finite for the leading order perturbative diagrams: 



F M (x,£,b A 



(1-x) 



N+l 



N M 

x]TQ 

3=1 



1-X 



■ (i - 06 x 



1-x 



,(1+0& J 



. (13) 

Below we shall explicitly compute this decomposition in 
the (J) 3 and Yukawa models — see Eqs. JSTJ, (J33JI, (pMjl . 

The triangle graph GPDs H M (x,^,t) obey the poly- 
nomiality constraints and automatically have the double 
distribution representation (|lip which naturally appears 
in terms of the a-parameter calculation of Feynman dia- 
grams HUG}. 

Now we can combine the physical and mathematical 
approaches. Triangle graphs will provide us with func- 
tions Q° M and with the corresponding double distribu- 
tions. Taking functions Q 3 M generated by the triangle 
graphs in various models, we can use these functions Q M 
in the general decomposition (|f Op . In this way we can 
construct GPDs obeying both polynomiality and posi- 
tivity constraints. 

The next step is to take perturbative theories contain- 
ing several parton fields with different masses. In this 
case asymmetric triangle graphs with different masses 
will enter the game. Taking the number of different 
masses to infinity one arrives at triangle graphs inte- 
grated over the masses. Under certain restrictions on 
the integration weight this will generate GPDs satisfying 
both positivity and polynomiality constraints. 

Another important ingredient is the D term JjjJ. For- 
mally one can use the trick of Ref. and include the D 
term in the double distribution representation for GPDs. 
However, for the analysis of the positivity bounds the 
explicit form of the D term is much more convenient. 
Indeed, the D term vanishes in the region \x\ > |£| and 
therefore it is not restricted by the positivity bound JSJ. 
On the other hand, the D term automatically satisfies 
the polynomiality constraint. This means that construct- 
ing the solutions of the polynomiality and positivity con- 
straints we are free to add an arbitrary D term. 
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FIG. 1: The triangle graph for GPDs. The light-cone mo- 
mentum fractions are normalized with respect to (Pi + P2)/2. 



V. TRIANGLE GRAPH IN THE 3 MODEL 

Let us start the analysis of the positivity properties 
of perturbative diagrams from the case of scalar partons. 
The perturbative triangle graph of the <f> 3 model of Fig. 
is often used as a toy model for GPDs Q ■ This triangle 
graph leads to the following Feynman integral: 



d 4 q 
(2tt)4 



2{nq) 



n{P 1 + P 2 



{Pi - q) 2 -m 2 + iOq 2 -m 2 + iO (P 2 - q) 2 - m 2 + iO ' 

(14) 

Here is the coupling constant of the interaction be- 
tween the scalar "partons" with the mass m q and the 
scalar meson with the mass M, 



Pi =P'i = M 2 . (15) 
It is assumed that the meson stability condition holds: 



M < 2m 



(16) 



We are going to consider the version of the (j> 3 model 
with such a flavor content and couplings which select the 
diagram of Fig. ^ but forbid the cross-channel diagram. 
The resulting GPD vanishes in the "antiparton region" 

*<-ie 

In Appendix[B]this diagram is computed in the general 
case of three different "parton" masses in the triangle. 
Setting these masses equal in Eq. I|B8|) 



(17) 



mi = "12 = n-3 = rn q . 



we obtain the double distribution representation for our 
case (|T4T) : 



«1 — ai 

H <t ,3(x,S,,t) = I dai / da 2 F < p{ai,a 2 ,t) 



where 



x<5[x-l + ai(l + e) + a 2 (l-0] 



(18) 



F^(ai,a2,*) 



9^ 

16tt 2 



x [m 2 - (ai + a 2 )(l - ai — a 2 )M 2 - aia 2 t] 



(19) 



The afc parameters used here are related to parameters 
a, (3 appearing in Eq. @ as follows: 

ai = i(l-/3-a), a 2 = l -(l-p + a). (20) 

In Appendix El the triangle graph is computed in the 
impact parameter representation in the region x > |£| for 
the case of three different "parton" masses. Setting these 
masses equal in Eq. (|C10ll we find the impact parameter 
representation © for our graph at x > |£| 



l-x 

47T 



x V 



V 



(21) 



with the function V given by Eq. i|C6(l in terms of the 
modified Bessel function Kq: 



V(r,c x ) 



g<t>z 

2nr 



r(l 



M 2 



(22) 

The factorizcd form of the result l|21|) for the GPD in 
the impact parameter representation obtained in the 4> 3 
model is an illustration of the general decomposition of 
triangle diagrams l|13|) . We see that in our case the sum 
on the RHS of Eq. (|13|) contains only one term. The rea- 
son is that the q propagator of our diagram corresponds 
to a spin-zero particle. 

Introducing the variables (see Appendix ^ for more 
details) 



1 



1 



n = 



r 2 = 



(23) 



instead of x, £, and working in the region x > |£| (i.e. 
< J"i,r 2 < 1), one can rewrite Eq. 121|) in the form 



r . {.r.a.-L^i, 



1 nr 2 

2n n + r 2 



V{ri,nb^)V{r 2 ,r 2 b^) {x > . 



(24) 
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FIG. 2: Reduced diagrams coming from the triangle graph in Yukawa model. 



VI. TRIANGLE GRAPH IN YUKAWA MODEL 



Now let us compute the "quark-in-meson" GPD in Yukawa model. The same triangle graph of Fig. (now with 
the fermion loop) leads to the following Feynman integral 



x — 1 



2(nq) 



n(Pi + P 2 ) 



x Tr 



(717) 



(Pi - 9)7 ~m q + iO -(qj) -nig + iO (P 2 - q)j 



iO 



(25) 



The factor of 1/2 on the RHS is inherited from the light-ray fermion operator O^ 1 (see TableHJ, and gy is the coupling 
constant. Again we assume that the flavor content and couplings are chosen so that the cross-channel triangle diagram 
is forbidden so that we deal with the GPD vanishing at x < — 1£|. 

The trace of the Dirac matrices can be represented in the following form: 



Tr {(717) [(Pi - q)j + m q ] [-(qj) + m q ] [(P 2 - q)j + m q }} 



-P 



7 Pa 



(nPi) [(P a - qf - ml] + rK 2 ) [(Pi - q) 2 - m*] -(nq) [Am 2 - (PP 2 )] - [n (Pi + P 2 )] 



-M 



2m 2 



(26) 



Most of the terms on the RHS of Eq. (|26|1 contain factors which cancel one of the propagators in the denominator 
so that one arrives at reduced diagrams containing only two propagators (Fig.(5J. The contribution of the nonreduced 
part is 



H Y {x,Z ) t)=2& J ^ 



x - 1 + 



2{nq) 



n(Pi + P 2 



l -(nq) [Am 2 - (PiP 2 )] - [n (P + P 2 )] Qm 2 - 2m 2 ) J 



(Pi - q) 2 - m 2 + iO q 2 - m 2 q + iO (P 2 - q) 2 - m 2 + iO 



reduced diagrams . 



(27) 



Comparing the RHS with Eq. (|14(l , we see that we have reduced the calculation of the GPD in the Yukawa model to 
the scalar GPD in the d> 3 model 



H Y (x,U) = 2 £ r- 



(1 - x)t - x (Am 2 - M 2 ) 



Hj,3 (x, £, t) + reduced diagrams . 



(28) 
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The reduced diagrams (a) and (b) of Fig. [5] give a t independent contribution. Therefore in the impact parameter 
representation they vanish if b ^ 0. The contribution of the reduced diagram (c) of Fig. [21 has the structure of a D 
term (JHJ) which vanishes ii \x\ > |£|. Thus all three reduced diagrams can be ignored if one is interested in the region 
\x\>\Z\,b L ?0. 

Let us transform Eq. (|28|l into the impact parameter representation omitting the reduced diagrams. Parameter t 
© becomes a differential operator: 



4£ 2 Af 2 



(— V 



4£ 2 M 2 



With this expression for t and with the representation (|21[1 for {x, £, b- 1 ) we find from Eq. 128|) 



2 1 - 



x (4to 2 - Af 2 ) 



(29) 



x V 



1 -a; 



>(i-0& J 



u 



i 



(30) 



Functions V (r, c x ) expressed in terms of the modified Bessel functions (12211 obey the following differential equation 

(a? 7 ) v ^ c± ^ = ^ rM2 + r_lm «) ~ m2 ] r ~ ly ( r > c± ) • ( 31 ) 

Using this differential equation and variables <|23[1 we find from Eq. (|3U[) 

" -xf4-ir^— H [(l-2r 1 )y(r 1 ,r 1 6 ± )] [(1 - 2r 2 )U(r 2l r 2 ^)] 



h- L #0,a;>|{| 



(32) 



We see that the RHS has the general structure 113(1 satisfying the positivity bounds. 

One can also compute the GPD ffy i76 (:r,£,i) in the Yukawa model with the pseudoscalar coupling. Replacing the 
interaction gy(f>ipip — * gy ^i^tp , one slightly changes the Dirac trace l|26|) . which leads to the following modification 
of the GPD: 



i9y 



HY, l5 (x,£,t) = H Y (x,£,t) +8m l} —xH < p3(x,£,t) 

9<f>3 



(33) 



In the impact parameter representation we again find an example of the general structure l|13[) which guarantees the 
positivity: 



1 — x 

F Y,tB ( x, f , 1 _ ^ 2 b A 



Fy ( x >Z>Y^2 bX ) +8m 2 x^F a ( x >Z>Y^2 b± 



(1 



^9 Y 



1% 27rrir 2 



{m 2 F(r 1 ,r 1 6 ± )7(r 2 ,r 2 6 ± ) + [V 6 xU (r Xj r x 6 x )] [V 6 x7 fa, fair 1 )]} 



r 



(34) 



VII. BUILDING MODELS OF GPDS FROM 
TRIANGLE GRAPHS 



In the previous section we have explicitly checked that 
the triangle diagrams in Yukawa model generate GPDs 



satisfying both polynomiality and positivity constraints. 
Since the fermion-in-(pseudo)scalar GPDs obey the same 
polynomiality and positivity constraints in the Yukawa 
model and in QCD we can use the triangle GPDs of the 
Yukawa model as elements for the construction of models 
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of the quark GPD in pion, compatible with the polyno- 
miality and positivity constraints. 

The first step is to mix the GPDs of the scalar and 
pseudoscalar Yukawa model 



C 1 H Y (x,^t)+C 2 H Yn5 (x,^t) 



(35) 



with positive coefficients, which is equivalent to the 
Yukawa model with the coupling 



(36) 



SyH> [y C\ + «75 v C2 ) tp ■ 
Now let us show that the function 

(1 - x)H^(x,^,t) 



(37) 



also obeys both polynomiality and positivity constraints 
for the quark GPD in pion. Indeed, the positivity in- 
equality JHl for the fermion- in-scalar GPD (N = 1) dif- 
fers from the case of the scalar-in-scalar GPD (N = 0) 
exactly by the factor of (1 — x). The polynomiality con- 
dition {JJ for the fermion-in-scalar GPD also allows one 
more degree of x compared to the GPD in the </> 3 model. 
Now we can use all available elements, l|35f) and l|37|l . 
to build models for the pion GPD. For any positive co- 
efficients Cfe the following combination will satisfy both 
polynomiality and positivity constraints: 



dH Y (x, £, t) + C 2 H Yn5 (x, f , t) + C 3 (l - x)H^ (a;, £, t) 



A K (Ci + C 2 ) 



-(1 - x)t - x (4m 2 q - M 2 ) 



9 2 ^ 



+4C 2 m 2 x + C 37 f T (l-x) \Hp(x,£,t). 



(38) 



The next step is to consider triangle graphs with arbitrary masses. We start from the </> 3 model. Let us take the 
triangle graph of Fig. ^ with the masses m\ for the (Pi — (j)-propagator, m 2 for [P% — q) and 771,3 for q. This graph is 
computed in Appendix [5] The result l|B8() can be represented in the following form: 



/>1 pl—OLX 

H lfj 3(x,£,,t\m 1 ,m 2 ,m 3 ) = / da>i / da 2 F® a (a 1 ,a 2 ,t\m 1 ,m 2 ,m 3 )6 [x - 1 + ^(1 + + a 2 (l - £)] , (39) 

Jo Jo 



FR(a 1 ,a 2 ,t\mi,m 2 ,m 3 ) 



167T 



{~aia 2 t - (ai + a 2 ) [1 - (ai + a 2 )} M 2 



\a\m\ + a 2 m 2 + (1 — ai — 0:2)7773] } 1 . 



(40) 



The corresponding impact parameter representation in the region x > |£| computed in Appendix [U] is given by 
Eq. iUTOl) : 



1 — x 
F^s [ x,£, - -6 J 



i-e 

The function V is defined by Eq. l(C6} 



my. iii 2 - in.- J = — (1 - x)V(ri,r%b \m\, m 3 )V(r 2 , r 2 b \m 2 ,m 3 ) . 



Vin^mums) = T^Ko [\^\r^ ^ ri m 2 + (1 - ri )m 2 - n (1 - n) M 2 ^j 



(41) 



(42) 



Any single triangle graph automatically satisfies the 
polynomiality constraint. Therefore mixing the contri- 
butions of various triangle graphs we must take care 
only about the positivity. Keeping in mind the factor- 
ized structure I|41J1 we see that the integral 

H {0 \x,£_,t) = J dmx J dm 2 J dm 3 



x s(mi,m2,m 3 )Hj > 3(x,£,t\mi,m 2 ,m 3 ) (43) 

is compatible with the positivity if the weight s has the 
structure 



s(mi, 777 2 , m 3 ) = / d\a(mi,m 3 ,\)a*(m 2 ,m 3 ,X) (44) 
where the function a(mi, 777,3, A) is arbitrary. This rep- 
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resentation is equivalent to the following property of 
s(m 1 ,m 2 ,m 3 ): 



dmi / dm 2 s(mi,m 2 ,m 3 )f(mi)f*(m 2 )>0 (45) 



for any function f(m) and for any value of 7713. Since we 
are interested in real GPDs even in £, we must work with 
real functions a&. 

It is also assumed that functions a(mi, 7773, A) are com- 
patible with the stability of the meson: 

a(mi,m3,A)=0 if mi + 7773 < M . (46) 



Now we can turn to Yukawa model, generalize the 
representation l(;-i8(l to the case of different masses 
?7Ji,m2, TO3 and integrate over these masses by analogy 
with the <fi 3 model, Eq. I|43l) . The generalization of 
Eqs. l(2"51) and for the case of different masses nik 
is 



H Y (x, £, i|77ii, 777 2 , m 3 ) = ( ) { (1 — x)t + 2xM 2 - 2x(m\ + m 3 )(m 2 + m 3 ) - (mi — 7772)^ 



+£(wi — 1712) (mi + m2 + 2m3)} Hj,3(x, £, t\mi, mi, m 3 ) + reduced diagrams , 



(47) 



HY, 15 (x,£,t\mi,m2,m3) = H Y (x,£,t) +4 



Combining the "superposition" 



9Y_ 

9<p 3 



m 3 [(x - £) m i + (x + Qmi] H <f ,3(x,^t\m 1 ,m 2l m 3 ) . (48) 



J dmi J dm 2 J dm 3 \ Oj^J [ s i( m i: m 2, m 3 )H Y (x, £,, t\mi , m 2 , m 3 ) + s 2 (m 1 , m 2 , m 3 )H Y ^ (x, £, t\m x , m 2 ,m 3 )] 



+ (1 - x)s2("7i,m 2 ,m 3 )i?^3(a;,^,t|mi,m 2 ,m3) 



(49) 



with a -D term and using Eqs. I|47|) . I|48|) . we arrive at the following solution of the positivity and polynomiality 
constraints for the fermion-in-scalar GPDs (i/W with N — 1): 

(#,£,*) = / dmi / dm 2 / dm 3 H^3 (x, £, t|mi, m 2 , m 3 ) {[si(mi, m 2 , m 3 ) + s 2 (mi, m 2 , m 3 )] 



x [(1 — x)t + 2a;M 2 — 1x{m\ + m 3 )(m 2 + m 3 ) —(mi — 7712) 2 + £(mi — m2)(mi + 2m 2 + 2m3)] 



+ 4s 2 (mi,m 2 ,m 3 )m3 [(x - £)mi + (x + £)m 2 ] +s 3 (mi, m 2 , m 3 )(l - x)} + 6 (|£| - |x|) D [ -,t) sign(£) . (50) 



r 



Here the integration weights must have the structure 



Sfe(mi,m 2 ,m 3 ) = / d\a k (m 1 ,m 3 ,\)a* k {m 2 ,m, 3 ,\) 



(51) 

with arbitrary functions a/ c (mi,m 3 ,A) obeying the me- 
son stability condition 1)46(1. The functions a k must be 
real if one is interested in real £-even GPDs. We remind 
the reader that the term D(x/£,, t) on the RHS of Eq. lj5U|) 



is not constrained by the polynomiality and positivity. 

The triangle GPD vanishes in the antiquark re- 
gion x < — 1£|. Therefore the construction i|50|l should 
be modified by adding a similar contribution with the 
replacement x — » — x and with its own set of coefficients 

Sk- 

As mentioned above, we have checked that the GPDs 
obtained from the triangle graphs satisfy the positivity 
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bounds in the impact parameter representation only at 
b ± ^ 0. At b 1 - = we must take into account the Sfi- 1 ) 
contributions coming from the reduced diagrams (a), (b) 
of Fig. [21 which depend on the normalization point /i and 
can violate the positivity bounds. This b = singularity 
of the triangle diagrams is the leading order perturbative 
manifestation of more serious problems which can be met 
due to a nontrivial interplay between the two scales 
and b 1 - EH . If one wants to construct models of GPDs 
avoiding this small b 1 - problem, then one can impose the 
following condition on the coefficients ak(Tn%, m 3 , A) ap- 
pearing in our construction of the integration weight 151|) 



<im 1 afe(m 1 , to 3 , A) = (A; = 1,2). 



(52) 



Indeed, the reduced diagram of Fig. |51 (b) does not de- 
pend on mi, therefore after the integration over the 
masses in Eq. H5()(l with the weight Q51JI obeying the con- 
dition Q52JI. the contribution of the diagram (b) vanishes. 
The contribution of the diagram (a) is 777,2 independent 
and vanishes after the integration over mi- Condition 
(|52|l also suppresses the unacceptable large t behavior of 
triangle graphs. Note that Eq. (|52|l means that the func- 
tions afe(mi, TO3, A) cannot be positive everywhere. This 



is not a problem because in order to satisfy the positivity 
bounds on GPDs we need only the construction 151fl for 
the functions Sk and we have no restrictions on the sign 
of the functions etfc. 



VIII. POSITIVITY OF FORWARD PARTON 
DISTRIBUTIONS 



The positivity of forward parton distributions (FPDs) 
is a consequence of the positivity bounds on GPDs. This 
idea is present in an explicit or implicit form practi- 
cally in all pap ers dea ling with the positivity bounds on 

gpds [n hi d q d e in d d m m mm 

. Since our construction of scalar-in-scalar GPDs (|43|) 
and fermion-in-scalar GPDs i|50|) satisfies the positivity 
bounds on GPDs (JHJ , the positivity of the corresponding 
FPDs is predetermined. Nevertheless the direct explicit 
check of the positivity of FPDs is rather interesting. In 
particular, in the case the fermion-in-scalar GPDs made 
of the triangle graphs of the Yukawa model, the analysis 
of the forward limit is instructive for understanding the 
role of the "divergence-cancellation" condition l|52|) . 



The FPD of the (f> 3 model is given by 



q <j> 3{x\m u m 2 ,m 3 ) = H <j) 3(x,^t\m 1 ,m 2l m 3 ) 



£=0, t=o 



9(x) 



I67H mf 



In 



xm\ 



+ (1 - x)m\ - x(l - x)M 2 



OTI3 + (1 — x)rri2 — x(l — x)M 2 



(53) 



This result can be obtained by taking the forward limit in Eq. i|39|) . Note that here we define the FPD q^z as the 
forward limit of the GPD H^z. This differs by a factor of x _1 from the physical definition in terms of the density 
of partons. If we take the FPD q^ for the (f> 3 model without integrating this FPD over masses mk, then in the 
symmetric case m\ = m 2 we find from Eq. (|53|l 



ps (x\m 1 ,m 1 ,m 3 ) 



9p 



(l-x)9{x) 



(54) 



(55) 



(56) 



16-7T 2 xm\ + (1 — x)m\ — x(l — x)M 2 
The positivity of this function is obvious if we assume the meson stability condition 

M < mi + 7713 ■ 

Indeed, 

xm\ + (1 - x)m\ - x(l - x)M 2 = [2:7773 - (1 - x)777i] 2 + x(l - x) [{m x + m 3 ) 2 - M 2 ] > . 

Next, we can consider the scalar- in-scalar GPDs constructed according to Eq. I|43|l and take the forward limit. With 
expressions l|44|) for s and l|53|) for q^ we find 

q(x) — J dmi J dm 2 J dm 3 s(mi 7 m 2 , m^q^s (x\mi, m 2 , m 3 ) 

6(x) j xm\ + (1 - x)m\ - x(l - x)M 2 



= I dX I dmi / dm2 / dm 3 a(mi, 777,3, \)a*{m 2l m 3 , A) 



167r 2 mf 



4 + (I- x)m\ 



x(l -x)M 2 



(57) 



The positivity of this expression is a consequence of the following inequality valid for any function / and for any v 
(see Appendix [D|) 



dmi 



dm 2 f{m 1 )f*(m 2 )- 



ln- 



> 0. 



HIT) 



(58) 



m\+v>0 
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Now let us turn to the Yukawa model. A straightforward calculation allows us to express the triangle graph 
contribution (with different parton masses) to the FPD of Yukawa model qy in terms of the triangle FPD q^i of the 
<fi 3 model as follows: 

q Y (x\m l7 m 2 ,m 3 ) = ( — ] \2xM 2 - 2x{m 1 + m 3 )(m 2 + m 3 ) - (mi - m 2 ) 2 ] qMx) 
\9<pJ 



g Y 0(x) ( (1 - x)m\ + xmj - x(l - x)M 2 



In 



(1 — x)m\ + xm\ — x(l — x)M 2 

7 2 



(59) 



The first term on the RHS containing q^ can be obtained by taking the forward limit in Eq. I|47|) . The logarithmic 
terms on the RHS are generated by the reduced diagrams (a), (b) of Fig. [3 The ultraviolet divergences of the Yukawa 
model are renormalized at the scale /x. For fixed parton masses, the FPD qy (|59|l depends on the normalization 
point (i via the additive term ln/x. This simple [i dependence obviously leads to the violation of the positivity at low 
normalization points and to the restoration of the positivity at large /x (here the formal behavior of the triangle graph 
is meant and not the properties of the full Yukawa model). 

Next we want to study the forward limit of the fermion- in-scalar GPDs constructed according to Eq. (|50l) . Since the 
logarithmic /x dependent terms in Eq. 1)59(1 depend either on mi or on m 2 but not on both rti\ and m 2 simultaneously, 
we conclude that these logarithmic terms will be cancelled by the integration over mi and m 2 due to the condition 
(|52p. For simplicity let us consider the case a 2 = a 3 = 0. Then Eq. I|5U|I generates the following FPD: 



dmi / dm 2 / dm 3 s 1 (mi,m 2 ,m 3 )q Y (x\rn 1 ,rn 2 ,m 3 ) 



167T 2 ' 



9(x) dX dmi / dm 2 I dm 3 ai(mi, m 3l A) 



X al(m 2 , m 3 , A) [2xM 2 - 2x{mi + m 3 )(m 2 + m 3 ) - {mi — m 2 y 



The positivity of this FPD reduces to the following inequality: 



In 



xm\ + (1 — x)m\ — x(l — x)M 2 
xm\ + (1 — x)m 2 — x(l — x)M 2 



(60) 



dX / dmi / dm 2 / dm 3 ai(mi, m 3 , X)a* 1 {m 2: m 3 , X)- 



x [2xM 2 — 2x(mi + m 3 )(m 2 + m 3 ) — (mi — TO2) 2 ] In 



xm\ + (1 — x)m\ — x(l — x)M 2 
xm\ + (1 — x)m 2 — x(l — x)M 2 



> 0. 



(61) 



Fixing A and m 3 , one can show that this inequality holds 
already after the integration over mi and m 2 . In order 
to see this, we first have to rearrange the factor in the 
brackets as follows: 



2xM 2 — 2x(mi + m 3 )(m 2 + m 3 ) — (mi — m 2 ) z 



Ai+A 2 



(62) 



where 
Ai = 2(1 - x) ( mi 



1 -a; 



m 3 m 2 



l-x 



m 3 



-4, 



njj + m 2 + 2x 



1 -x 



M* 



(63) 
(64) 



The positivity of the Ai contribution to the inequal- 
ity l|61|l follows from the inequality l|58|l combined with 



the meson stability condition (|46|l for the function 01. In 
order to prove the positivity of the contribution of A 2 to 
the inequality l|61|l . one has to use the following inequal- 
ity derived in Appendix iDl 



dmif(mi) 



dm 2 f*(m 2 ) 



2v m\ 



■hi 



<0. 



(65) 



This inequality holds for any function / obeying the con- 
dition 



J dmf(m) = . 



(66) 



In the case of the A 2 contribution to the inequality 161|l , 
the condition l|66|) holds due to Eq. I|52|) . 
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IX. CONCLUSIONS 



This work was supported by DFG and BMBF. 



In this paper, it is shown that the representation l|50(l 
for the quark-in-pion GPDs automatically obeys both 
polynomiality and positivity constraints. This construc- 
tion is based on the integration of the triangle graphs for 
Yukawa model over the masses of the three propagators. 
It also contains the piece (1 — x)H^3 whose positivity 
and polynomiality properties are inherited from the tri- 
angle graph of the 3 model. The mass integration al- 
lows a wide class of mass dependent weights constrained 
only by the positivity condition l|51[) . by the divergence- 
cancellation requirement (|52() and by the meson stability 
condition (|46|) . 

We also have the freedom of adding an arbitrary D 
term without violating positivity and polynomiality. The 
possibility to include the D term is very important. In- 
deed, integrating over the masses of triangle graphs one 
can generate only thresholds in the t channel whereas the 
D term allows us to produce single-particle poles in the 
t channel. 

This paper describes only the method of the construc- 
tion of GPDs obeying polynomiality and positivity con- 
straints. One can go beyond the </> 3 and Yukawa models 
trying to find new "perturbative bricks" for the construc- 
tion of the solutions of the positivity and polynomial- 
ity constraints. One should keep in mind that in con- 
trast to the two-point correlation functions for which we 
have Kallen-Lehmann representation, the case of GPDs 
is more involved and there is no guarantee that the true 
physical GPD can be represented as an integral of tri- 
angle graphs over their masses even if we go beyond the 
Yukawa model, include triangle graphs from other the- 
ories and make our best from the freedom to add an 
arbitrary D term. 

On the other hand, our construction is parametrized 
by arbitrary [up to the constraints (|46|l . l|52|l ] functions 
afe(mi,TO3,A) depending on three variables, i.e. our 
parametrization has the same amount of "degrees of free- 
dom" as the GPD H (x, £, t) which also depends on three 
variables. This means that the set of the solutions of the 
positivity and polynomiality constraints covered by the 
representation H5()(l is rather large. 

The comparison of the triangle graph approach con- 
sidered here with the formal mathematical solution of 
the positivity and polynomiality constraints suggested in 
Ref. [24j shows a number of similar features but at the 
moment it is not clear how large the overlap between the 
two representations is. As long as this issue is not clar- 
ified it makes sense to work with the "superposition" of 
the two representations. Indeed, from the practical point 
of view the variety of structures compatible with the 
polynomiality and positivity is more important than the 
problem of the unambiguous parametrization of GPDs. 

Acknowledgments. I am grateful to Ya.I. Azimov, 
A.V. Belitsky, M. Dichl, L. Frankfurt, D.S. Hwang, 
R. Jakob, P. Kroll, D. Miiller, M.V. Polyakov, 
A.V. Radyushkin and M. Strikman for useful discussions. 



APPENDIX A: SOLUTION OF THE POSITIVITY 
BOUNDS 

In this appendix we describe the properties of the vari- 
ables r\ , r 2 and derive the solution l|l()(l of the positivity 
bound ©• 

The variables r%, r 2 , which can be used instead of x, £, 
arc defined as follows: 



1 



'I = 



1 + 



1-C 



r 2 - ri 

T2 + J*l 



2rir 2 
r\ + r 2 



2dxd£ dridr 2 



(Al) 



(A2) 



(A3) 



(1 — x) 3 r\r\ 

The region covered by the positivity bound JSJl 

x > |C| (A4) 

is mapped onto the square in the r\ , r 2 plane 

<n,r 2 < 1. (A5) 

Inequality ijSJ) takes the following form in terms of inte- 
gration variables ri,r 2 (we keep notation £,£ in GPDs 
implying that these variables are functions of r 1; r 2 ) 



1 dri f 1 dr 2 f n+r 2 \ N+1 . , . 

-JT / -X — — J P (ra)p(ri) 

r l JO r 2 V r l r 2 / 



1 - x 



b ± \>0 



(A6) 



Since function p is arbitrary we can replace it 

p(n) - rf +3 P ( ri ) , (A7) 
which leads us to the equivalent form of inequality l|A6|l 



i pi 

dn dr 2 {ri+r 2 ) N+1 p* {r 2 )p(r{) 
o Jo 



1 - X 



> 



Inequality (|A6(I means that the function 

N + l 



n + r 2 
r\r 2 



(A8) 



(A9) 
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must be a positive definite quadratic form, i.e. it has the 
following representation 

(^f V> (*^) 

= Y / Rn(n,b ± )R* n (r 2 ,b ± ) (A10) 

n 

with some functions R n . Turning back to the variables 
x, £ we find 

(x,£,b ± )=(l-x) N+1 



APPENDIX B: DOUBLE DISTRIBUTION IN 
THE 3 MODEL 

This appendix contains a brief derivation of the double 
distribution representation 0,0)0] for the GPD in the 4> 3 
model. The contribution of the triangle graph of Fig. 
with three different masses of partons is 

H <t> 3(x, £,t\mi,rri2,ma) = {ig^) 2 

X / 1 TT^ \X - H x 

y (2tt)4 v P+ + P+J 



(All) 

In the case of real and £-even GPDs, the functions R n 
are real. Introducing the functions 

Q n (r,b^)=R n (r,h ± ^ , (A12) 
we obtain representation i|l(J|l for pW (x, £, 6^). 



i i i 

(Pi - q) 2 -mf + iO (P 2 - q) 2 - m\ + iO q 2 - m§ + iO ' 

(Bl) 

The light-cone components J 4 ± of vectors A 11 are assumed 
to be chosen so that the vector n appearing in the defini- 
tion of GPDs Q has only one nonvanishing component 

rT ": 

n+ = , n x = . (B2) 
Using the standard Feynman trick 



3 • pi p 1 — a i /• oo [ j 

TT — § — = / / ^ a 2 / d\ A 2 exp < iAV^ o*-k{q\ — m 2 + iO) > (B3) 

^—i \ h, — i / a 3 -l-a 1 -a 2 

with 

q 1 =P 1 -q, q 2 = P 2 -q, 93 = 9, P? = P§ = M 2 , (B4) 

we find 

rl rl-cti poo p / 2q^ 

H 4 ,3(x,£,,t\m 1 ,m2,m 3 ) — {ig^) 2 / dai / da 2 / d\X 2 4 <5 I x - 1 + 



o Jo 



(2tt) 4 V A + +^2 



x exp (iA {ai [(Pi - qf - m\\ + a 2 [(P 2 - ?) 2 - m 2 ] + (1 - ai - a 2 ) (g 2 - m 2 ) }) . (B5) 
The calculation of the integrals over q and A is straightforward and yields 

r dx a2 / i x ~ 1 + p^+w) exp (a {ai [(pi ~ 9)2 - m?] + aa [{P2 qr ~ m ^ 

+ (1 - ai - a 2 ) (, 2 - m 2 )}) = * 5 f * - 1 + 2°^ + ^ 



16tt 2 V p+ + P+ 

x {aia 2 (Pi - P 2 ) 2 + («i + a 2 ) [1 - (ai + a 2 )] Af 2 - [aim 2 + a^m 2 . + (1 - a x - 0^2)^3] } 1 • (B6) 
Taking into account that according to Eqs. Q, l|B2(l 
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we find from Eqs. (|B5f) . (|B6f) 

-ff03(x,^,i|mi,m 2 ,TO 3 ) 



503 



1 pi — OL\ 

da t I da 2 5 [x - 1 + ai(l + £) + a 2 (l - £)] 
o Jo 



x { [a^m 2 + a^m 2 + (1 — Q!i — o^to 2 ] — («i + a 2 ) [1 — (ai + 02)] M 2 — a^a^t} 



(B8) 



APPENDIX C: IMPACT PARAMETER 
REPRESENTATION FOR THE GPD IN THE 3 
MODEL 

In this appendix we compute the triangle graph for the 
GPD of the 4> 3 model in the impact parameter represen- 
tation in the region x > |£|. In principle this could be 
done by applying the Fourier transformation J^J) to the 
double distribution representation for this GPD dB8|) . 

However, we prefer another method based on the di- 
rect calculation of the diagram in the impact parameter 
representation (see e.g. Ref. The advantage of this 



r 



approach is that it explains the origin of the factorized 
form of the result. 

We start from the Feynman integral IjBljl for the trian- 
gle graph of Fig.^ One can integrate over q~ deforming 
the integration contour. At x > |£| the integral is deter- 
mined by the residue of the pole at q 2 = m 2 . Therefore 



at x > |£| we can replace on the RHS of Eq. (|BT 



2^{q+)5{q 2 - to 2 ) . (CI) 



q 2 — mg + iO 



Then we have 



H c j > 3{x,i,t\m 1 ,m,2,m i ) = (ig^f 



d 4 q 
(2tt) 4 



S[x-1+ 



P+ + P+ 



2tt% + )% 2 - TO 



2„_L 



d 2 q 



{Pi - q) 2 - to 2 + iO (P 2 - q) 2 -m 2 + iO 2(1 - x) J (2tt) 3 (Pi - q) 2 - to 2 + iO (P 2 - q) 2 - to 2 + iO' 
On the RHS the components q + and q~ are determined by the following equations 



Q = 



K + Pi 



(1-x), q 2 =m\. 



It is straightforward to show that 

(P k - q) 2 - to 2 = M 2 + to 2 - to 2 - 2(P k q) = M 2 + to 2 - to 2 - r k \p£ - ^7 Vf - (M 2 r k + m\rl 
Parameters r k are defined by Eq. I|A1(I . Using Eq. l|C4jl . we can rewrite Eq. I|C2|) as follows: 



H <j) 3(x,^,t\mi,m2,m 3 ) 



d 2 g A 

2{l-x) ri r 2 J (2tt) 



9^ 



(C2) 



(C3) 



(C4) 



IPj 1 -rvVI +TO^rf 1 (rf 1 - 1) + m 2 ^ 1 - M 2 (rf 1 - 1) 



Now wc define 



|P a x - r 2 ~V| + to 2 ^- 1 - ljrj 1 + m 2 ^ 1 - M 2 ^ 1 - 1) 



94> 3 f d2 V L , p ^b ± 



V {r k ,b m k ,m 3 ) = \ — 
r k J 2tt 



) 2 \ P ^\ 2 + mlrfir- 1 - 1) + to 2 ^ 1 - M 2 (rr 1 - 1) 



2nr k 



-K, 



(V kfc 1 ^rn 2 (\ - r k ) + m 2 r k - M 2 r k (l - r k )^j 



(C5) 



(C6) 
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where Kq is the modified Bessel function. Then it follows from Eq. I|C5(1 



■2 „2 



H rj>3 (x,£,t\m 1 ,m 2 ,m 3 ) 



47r(l - x) 

In the frame where P± + P 2 = 0, we have 



(fib 1 - exp [ib 1 - (r 1 P 1 L - r 2 P 2 )~\ V (r 1 ,r 1 b- L \m 1 ,m 3 ) V (r 2 , r 2 6 J -|m 2 , m 3 ) 



(C7) 
(C8) 



so that 

(x,£,t\mi,m 2 ,m 3 ) 



47r(l - a;) 



d 2 b exp 



V (ri,ri6 J -|mi,m 3 ) F (r 2 , r 2 6 ± |m2, m 3 ) . (C9) 



We see that in the impact parameter representation (jSJ our triangle graph for the GPD of the tfi 3 model has the 
following form: 



/ 1 — x 

lx,£, _ -$ b | mi, 



1 - .T 



m 2 ,m 3 ) = V" (ri,ri& |rai,m 3 ) V (r 2 ,r 2 b^\m 2 ,m 3 ) (x>|£|). (CIO) 



APPENDIX D: USEFUL INEQUALITIES 



In this appendix we derive two inequalities used in Sec. IVIIII 
Inequality 1. For any function / and for any real constant v 



dm x \ dm 2 f(m 1 )f*(m 2 ) 



1 , tot + v 
—5 o In —g > 

mf — to. Too + v 



Proof. Obviously 



dmi / c?to 2 /(toi)/*(to 2 ) 



m\ + v 

TTl? — TOo TOo + v 



hi 



dmi 



dm 2 f(m 1 )f*(m 2 ) 



mj+i>>0 m|+y>0 



(Dl) 



' d7 (TO?+7 + ^)(TO.|+7 + ^) ' 



mJ+i/>0 



dmi f (mi) 2 j ; — 

TO| + 7 + V 



> 0. 



Inequality 2. For any function /, obeying the condition 

dmi f (mi) = 



with some real we have 



dmi f (mi) J dm 2 f*(m 2 ) 

mi+v>0 



to? + TOo + 2v , to? + z/ 
In- 



< 0. 



Proof. Omitting for brevity the integration region m k + ;/ > 0, we can write using Eq. |jD 
dmi f [mi ) I dm 2 f (m 2 ) — g 



(D2) 



(D3) 



(D4) 



to? — m; 



. m-i + v f , . f , N rni + u mi + v 

m — ^ = zKe / ami j [mi) / dm 2 j (m 2 ) — 5 5 in — 5 

to 2 + v 7 7 to^ — m 2 to 2 + v 



= 2I{e J o d l J 



dmi 



m\ + v 
7 + m\ + v 



f{mi) 



dm 2 — ■ — — /* (to 2 ) 

7 + TO-2 + V 
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2Re / d-y 



dmif(mi) + / dm 



ra\ + v 
' 7 + m\ + v 



/(mi) 



dm 2 — ; — t — f*(m 2 ) 

7 + JTT-2 + V 



-2Re / d-f 



dm i 



7 + 772^ + V 



dmo 



f*(m 2 ) 
' 7 + m| + 2/ 



-2 / c?77 
'o 



G?TOl 



/(mi) 



' 7 + TO j + V 



< 0. 



(D5) 
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